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ABSTRACT 

The  objective  of  the  work  is  to  identify  the  flow  mechanisms  that  cause  the  decrease  in 
spreading  rate  of  supersonic  mixing  layers  as  the  convective  Mach  number  increases  and  to 
suggest  means  of  enhancing  the  mixing.  Two  approaches  have  been  taken,  one  numerical 
and  one  analytic.  A  computer  code,  TMRC,  has  been  used  to  simulate  both  time  and  space 
developing  mixing  layers  to  get  some  indication  of  the  flow  physics.  To  complement  the 
numerical  study  a  simple  analysis  has  been  developed  which  explains  the  variation  of 
mixing  rate  with  convective  Mach  number.  The  analysis  seems  to  indicate  that  little  can  be 
done  to  enhance  mixing  as  such  although  the  real  problem  of  simultaneous  mixing  and 
combustion  may  be  more  amenable  to  control. 

Research  Objectives 

The  overall  objective  is  to  identify  the  flow  mechanisms  that  cause  the  decrease  in  the 
spreading  rate  of  supersonic  mixing  layers  and  to  suggest  means  of  enhancing  the  mixing. 

In  the  past  year  the  objective  hajnrbeen  to  extend  and  consolidate  the  simple  theory 
developed  in  the  first  year  of  the  contract  and  to  see  if  it  provides  some  insight  in 
enhancement  of  mixing.  This  study  is  augmented  by  the  development  of  a  computer  code 
to  numerically  simulate  the  temporary  and  spatially  developing  mixing  layers  to  gain 
insight  into  the  physics  of  the  flow. 

Status  of  the  Research  Effort 

At  the  end  of  the  first  year  of  the  contract  a  simple  theory  had  been  developed  of  the 
spreading  rate  variation  with  convective  Mach  number.  This  theory  is  based  on  the 
transonic  small  disturbance  equation  (Ref.  1) 

(^2  -  +  0yy  =  M^  +  2M^  (f>^^  (1) 

where  0  is  a  perturbation  velocity  potential,  t  is  time  nondimensionalized  by  c/U^  where  c 
and  Uj,  are  characteristic  length  and  convective  velocity  respectively,  and  M^  is  the 
convective  Mach  number. 
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k  =  (7+1  )M^ 

The  nondimensionalized  mixing,  m,  is  related  to  its  value  at  incompressible  (M^,  =  0) 
mixing,  mj,  by 
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where  m^  is  the  nondimensional  mixing  associated  with  a  rotational  kinetic  energy  (RKE) 
of  (RKE)^  at  a  convective  Mach  number  of  M^. 

If  «  1,  the  Prandtl-Glauert  equation  (k  =  0)  may  be  used  and  Equation  (3)  becomes 


—  =  p^/(l  -  M^/2)  (4) 
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The  results  predicted  by  this  theory  show  good  agreement  with  experimental  data,  see 
Figure  1,  especially  for  <  0.7. 

At  the  end  of  the  first  reporting  period  it  was  concluded  that  mixing  is  reduced  because 
the  eddies  produced  a  decrease  in  size  as  increases.  In  the  last  year  it  has  been  shown 
that  this  conclusion  is  wrong. 

Over  the  last  year  the  theory  has  been  consolidated,  that  is,  fewer  assumptions  are 
made,  and  has  been  extended  to  include  three-dimensional  effects.  A  full  description  of  the 
theory  is  given  in  the  Appendix,  which  is  AIAA  Paper  90-0706.  A  summary  of  the  main 
findings  is  given  below. 


(a)  No  matter  the  size  or  number  of  spanwise  vortices,  the  mixing  will  decrease  with 
increasing  at  the  same  rate.  This  indicates  that  two-dimensional  eddy  break  up 
devices  may  not  work. 


(b) 


Streamwise  vortices  will  not  be  affected  by  to  a  first  approximation  and  the  mixing 
associated  with  these  vortices  is  unchanged  with  M^,. 


(c)  If  the  hypothesis  is  made  that  large  fluid  mechanical  structures,  such  as  vortices,  will 
deform  to  avoid  shock  waves,  then  the  theory  is  greatly  improved  (Figure  1).  The 
deformation  may  consist  of  vortex  sweep,  rather  like  a  swept  wing  or  a  bifurcation, 
rather  like  the  braids  commonly  seen  in  incompressible  mixing  layers.  This  hypothesis, 
if  true,  would  have  considerable  implications  in  other  fluid  mechanics  problems. 
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(d)  Since  the  present  theory  is  based  on  notions  of  irrotational  flow  it  suggests  that  other 
studies  based  on  vorticity  redistribution  arguments  are  erroneous.  It  may  also  be 
suggested  that  since  the  present  theory  is  based  on  nonlinear  finite  disturbances, 
explanations  of  the  phenomena  based  on  linear  stability  theory  are  irrelevant.  This 
latter  viewpoint  is  corroborated  by  the  fact  that  experimental  data  are  taken  a 
considerable  distance  downstream  of  the  region  where  linear  instability  theories  are 
valid. 

(e)  The  main  cause  of  the  decrease  in  mixing  as  M^,  increases  is  that  an  increasing  amount 
of  the  energy  in  the  two  streams  goes  into  compression  of  the  fluid  during  entrainment 
rather  than  being  used  to  produce  pure  mixing.  The  only  way  to  increase  mixing  is  to 
accentuate  streamwise  vorticity  and,  if  the  hypothesis  in  (c)  is  valid.  Nature  is  doing  its 
best  to  achieve  this  goal. 

In  addition  to  the  theory  given  in  the  Appendix  numerical  simulations  of  both  time  and 
space  evolving  mixing  layers  have  been  performed.  The  code  used  is  TMRC  which  is  a 
fourth  order  accurate  code  developed  at  Nielsen  Engineering  &  Research  (NEAR)  for 
impinging  jet  flows.  These  calculations  were  performed  to  get  some  insight  into  two- 
dimensional  mixing  layers;  a  typical  example  of  such  a  calculation  is  shown  in  Figure  2.  At 
present  the  numerical  results  are  similar  to  those  obtained  by  other  investigators.  Attempts 
to  avoid  duplication  by  obtaining  the  computer  code  developed  by  NASA  Ames  Research 
Center  to  study  mixing  layers  proved  fruitless  on  the  grounds  that  the  code  was  not  ready 
for  public  use. 

Written  Publications 

Nixon,  D.,  Keefe,  L.  R.,  and  Kuhn,  G.  D.,  The  Effects  of  Compressibility  on  a  Supersonic 
Mixing  Layer,  AIAA  Paper  90-0706, 1990. 

Professional  Personnel 

Dr.  David  Nixon 
Dr.  Laurence  R.  Keefe 
Dr.  Gary  D.  Kuhn 

Interactions 

(a)  An  early  version  of  the  theory  was  presented  at  the  AFOSR/ONR  contractors  meeting 
in  Ann  Arbor,  Michigan  in  June  1989. 

(b)  Extensive  discussions  about  the  theory  were  conducted  with  Professor  Peter  Bradshaw 
of  Stanford  University,  June  1989. 

(c)  Discussions  on  the  theory  with  Dr.  Richard  Rivir  of  Air  Force  Propulsion  Laboratory, 
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(d)  Invited  seminar  on  theory  to  Minnesota  Supercomputer  Center,  October  18, 1989; 
cancelled  on  day  of  lecture  due  to  San  Francisco  earthquake;  rescheduled  for  April  25, 
1990. 

(e)  Paper  on  theory  to  be  presented  at  AIAA  Aerospace  Sciences  Meeting,  Reno,  Nevada, 
January  1990. 

(f)  Discussion  on  validity  of  the  hypothesis  noted  in  section  (c)  of  "Status  of  the  Research 
Effort"  with  Professor  Max  Platzer  of  Naval  Postgraduate  School,  Monterey,  CA,  with 
emphasis  on  relevance  to  dynamic  stall,  November  1989. 

(g)  Discussions  in  October  1989  on  validity  of  hypothesis  in  section  (c)  of  "Status  of  the 
Research  Effort"  and  other  aspects  of  mixing  layers  with  Mr.  Richard  Dix  of  AEDC, 
Tullahoma,  TN,  with  application  to  the  flow  in  weapons  bays. 

Discoveries 

No  new  discoveries  or  patents  have  resulted  from  this  work. 

Reference 

1.  Nixon,  D.  (editor):  Transonic  Aerodynamics,  Prog,  in  Astro  and  Aero.,  Vol.  81, 1983. 

2.  Papamoschou,  D.  and  Roshko,  A.;  Observation  of  Supersonic  Free  Shear  Layers,  AIAA 
86-0162,  January  1986,  Reno,  Nevada. 


-4- 


Figure  1.  Experimentally  measured  layer  growth  rates  compared  to  two-  and  three- 
dimensional  theory. 


Figure  2.  Vorticity  contours  in  a  two-dimensional  mixing  layer.  Mi  =  2.0,  M2  =  1.2, 
Me  =  .4  . 
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ABSTRACT 

A  theory  that  predicts  the  observed  decrease  of 
mixing  in  a  compressible  shear  layer  with  respect  to  its 
incompressible  counterpart  at  the  same  density  and 
velocity  ratios  is  described.  Unlike  previous 
investigations,  it  does  not  attribute  this  decrease  to  the 
reduction  in  growth  rates  calculated  from  linear 
instability  analysis  of  compressible  shear  layers.  Rather 
a  two  part  transonic  theory  is  offered.  At  lower 
convective  Mach  numbers,  the  asymptotic  mixing  rate  is 
calculated  analytically  by  relating  it  to  the  mixing  and 
energy  of  an  isolated  vortex,  and  then  replacing  that 
with  the  bound  circulation  of  a  two-dimensional 
transonic  airfoil.  At  higher  convective  Mach  numbers 
the  mixing  rate  is  well  described  by  assuming  an 
increasing  three-dimensionality  of  the  vortices  that 
delays  the  appearance  of  shock  waves.  This  three- 
dimensionality  can  be  represented  either  by  swept 
vortices  or  by  assumption  of  a  bifurcation  in  which 
longitudinal  vortices  become  increasingly  important. 

INTRODUCTION 

Compressible  mixing  layers  barely  live  up  to  their 
name.  To  the  dismay  of  supersonic  combustor 
designers,  adjacent  streams  of  differing,  high  speed, 
compressible  fluids  mix  much  slower  at  their  interface 
than  the  same  fluids  at  low  speeds  and  the  same  velocity 
and  density  ratios.  This  effect  was  demonstrated 
conclusively  in  the  experiments  of  Papamoschou  and 
Roshko,’  who  correlated  shear  layer  spreading  rate  data 
from  theirs  and  earlier  studies  using  the  convective 
Mach  number  (M^. )  concept  introduced  by  Bogdanoff.  ^ 
Since  then  other  experiments^  ,4,5  further 

documented  this  effect  both  in  plane  and  circular  shear 
layers.  Density  differences  between  the  streams  were 
originally  suggested  as  the  cause  of  this  mixing  decrease, 
but  the  work  of  Brown  and  Roshko*  showed  this  not  to 
be  the  case.  More  recently  this  phenomenon  has  been 
attributed  to  the  decrease  with  Mach  number  of  linear 
instability  growth  rates  calculated  for  supersonic  shear 
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layers.’ Lele®  has  provided  a  physical  explanation  for 
the  attenuation  of  linear  instability  by  pointing  out  the 
retarding  effects  of  velocity  divergence  and  baroclinic 
torques  in  the  initial  vorticity  redistribution  process. 

While  the  similarities  between  linear  growth 

reduction  and  mixing  rate  decrease  are  suggestive, 

positing  a  cause  and  effect  relation  between  them  ignores 

a  simple  experimental  fact;  measured  mixing  rates  are 

obtained  in  fully  turbulent  regions  of  the  flow,  far 

downstream  of  stations  where  linear  instability  is  the 

dominant  mechanism.  The  often  cited  measurements  of 

Papamoschou  and  Roshko  derived  mixing  rates  from 

stations  greater  than  1000  6,  where  6  is  the  trailing  edge 

boundary  layer  thickness  on  the  high  speed  side.  For 

this  region  they  state  (p.  462,  Ref.  1):  "Given  that  the 

shear  layer  is  fully  turbulent,  we  expect  it  to  grow 

linearly."  The  recent  measurements  of  Samimy  and 

Elliot^  were  obtained  from  flows  where  the  trailing  edge 

boundary  layers  were  fully  turbulent,  and  still  they 

deduced  the  mixing  rates  from  stations  at  least  twice  the 

distance  downstream  from  where  their  measurements 

show  the  turbulent  shear  layer  to  be  fully  established. 

The  spatially  developing  numerical  simulations  of  Lele 

begin  with  a  forced  laminar  tank  profile,  whose 

development  is  only  followed  for  some  200  initial  layer 

thicknesses  6  .  Growth  rates  are  display- 

w 

ed  for  only  the  initial  60  6  before  pair- 

u 

ing  even  occurs.  Thus  while  instability  theory  and 
experiments  show  similar  trends,  the  flow  regimes  of  the 
two  cases  are  not  the  same. 

In  the  following  sections  we  describe  a  theory 
appropriate  for  the  asymptotic  mixing  rate  of 
compressible  shear  layers;  it  does  not  depend  on  linear 
instability  arguments.  We  first  derive  an  expression  for 
the  energy  within  a  control  volume  circumscribing  a 
portion  of  the  turbulent  shear  layer,  in  terms  of  the 
mixing  rate  (rate  of  increase  of  mass  within  the  volume) 
and  convective  Mach  number  M^ .  We  then  show  that 
the  mixing  rate  of  an  isolated  vortex  is  proportional  to 
the  square  root  of  its  energy  and  thus  establish  a  scaling 
between  the  mixing  rate  due  to  an  arbitrary  vortex  and 
that  of  a  reference  vortex,  based  on  the  ratio  of  their 
energies.  Identifying  the  energy  and  mixing  rate  of  the 


control  volume  analysis  with  that  of  an  arbitrary  vortex 
yields  a  relation  between  the  nondimensionalized 
mixing  in  the  shear  layer  and  the  energy  and  mixing  due 
to  a  reference,  isolated  vortex  sitting  in  a  uniform  free 
stream  of  Mach  number  M^. 

The  reference  vortex  is  then  taken  to  be  the  bound 
circulation  of  a  lifting  airfoil.  Both  the  mixing  and 
energy  of  this  reference  "vortex"  can  be  calculated 
analytically  from  solution  of  the  Prandtl-Glauert 
equation  with  =  M^. .  The  resulting  expression  for 
the  shear  layer  mixing  as  a  function  of  agrees  well 
with  experimental  data  for  M^.  s  M*  -  .77.  Above  this 
critical  value  of  M^.  we  theorize  that  three-dimensional 
effects  become  more  important,  and  that  the  flow 
evolves  to  avoid  shocks.  One  configuration  that 
accomplishes  this  assumes  the  spanwise  vortices  acquire 
a  sweep  angle  just  sufficient  to  keep  velocities  subsonic. 
Another  possible  configuration  echoes  the  roller  and 
braid  structure  of  incompressible  layers.  Here  the 
spanwise  vortex  decays  to  avoid  shocks,  but  the 
streamwise  vortex  increases  in  strength.  In  either  case, 
simple  modifications  to  the  original  analysis  show  the 
normalized  mixing  rate  to  remain  essentially  constant  at 
its  value  at  M’,  again  in  agreement  with  experimental 
data. 


The  physical  configuration  of  a  compressible  free 
shear  layer  is  shown  in  Figure  1,  where  vorticity 
contours  have  been  drawn  from  the  results  of  our  own 
two-dimensional  direct  numerical  simulations  of  this 
flow.  These  simulations  integrate  the  Navier-Stokes 
equations  using  standard  fourth-order  finite  volume 
spatial  discretization,  fourth-order  Runge-Kutta  time 
advance,  and  "characteristic"  boundary  conditions  due  to 
Thompson^  °  at  the  outflow  and  upper  and  lower 
boundaries.  Two  parallel  streams,  initially  separated  by 
a  splitter  plate,  and  semi-infinite  in  directions  transverse 
to  the  mean  flow,  mix  at  their  interface.  The  upper 
stream  is  characterized  by  pj ,  U.| ,  a^ ,  and  the  lower  by 
Pj ,  Uj ,  aj  (Uj  2  Uj  assumed).  Static  pressures  in  the 
stream  are  matched  at  the  trailing  edge  of  the  splitter 
plate  so  that  no  expansion/compression  waves 
contaminate  the  streamwise  development  of  the  flow. 
Initial  instability  of  this  shear  layer  causes  it  to  roll  up 
into  spanwise  vortices  which  subsequently  pair  due  to  a 
further  spatial  subharmonic  instability.  As  the  roll-up 
and  growth  of  vortices  proceeds,  the  width  of  the  shear 
layer  grows  downstream,  and  increasing  amounts  of 
fluid  are  entrained  from  both  sides.  The  growth  rate  of 
incompressible  shear  layers  depends  on  their  velocity 
and  density  ratios.  High  speed  shear  layers,  however, 
grow  much  more  slowly  than  their  incompressible 
counterparts  at  the  same  velocity  and  density  ratios. 
Bogdanoff,  and  later  Papamoschou  and  Roshko, 
suggested  that  the  effect  of  compressibility  on  growth  of 
a  shear  layer  is  characterized  by  Mach  num¬ 
bers  M  measured  relative  to  the  convec- 
c , 

tivc  speed  tl^  of  the  large  structures.  Thus 


M 


-  u 
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M 

=2 


U 
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-  U 
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(1) 


A  relation  between  M  and  M  is  ob- 

tained  by  assuming  that  the  two  streams  have  equal 
static  pressures,  and  that  their  stagnation  values  in  the 
reference  frame  are  also  equal.  If  it  is  further 
assumed  that  the  ratio  of  specific  heats  7;  for  each 
stream  is  equal,  then  M^*=  “  M^,  and  the 

convective  velocity  may  be  calculated  as  a  sound  sfjeed 
weighted  average 


U 

c 


^2^1  ^1^2 

a  +  a 
1  2 


Although  there  is  some  scatter,  all  experimental  evidence 
supports  this  claim,  and  numerical  simulations  by  Lele  ^ 
have  shown  that  pressure  maxima  (stagnation  points) 
and  minima  (structure  centers)  do  travel  at  in  the 
initial  part  of  the  layer.  The  decrease  in  growth  rate  as 
Me  increases  can  be  seen  in  Figure  2,  where  the  thickness 
of  three  shear  layers  at  =  .4,  .6,  .7  are  plotted  versus 
downstream  distance.  The  width  measure  used  is  the 
vorticity  thickness. 
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CONTROL  VOLUME  ANALYSIS 


Refer  again  to  Figure  1  and  the  control  volume 
drawn  about  a  portion  of  the  flow.  The  volume  moves  at 
U^. ,  the  velocity  of  the  large  structures,  and  is  assumed  to 
contain  at  least  one  of  them.  The  upper  and  lower 
boundaries  are  sufficiently  far  from  the  interface  that  we 
assume  no  vertical  velocity  on  them.  Then  the  idealized 
model  of  the  mixing  layer  is  shown  in  Figure  3,  where 
flow  variables  U,  p,  p,  etc.,  are  subscripted  according  to 
the  volume  (Vj  ,  Vj  )  they  are  associated  with,  and 
quantities  topped  by  a  tilde  (-)  refer  to  average 
quantities  in  similarity  variables  on  the  segments  AD 
and  BC. 

In  Figure  3  the  difference  between  fluid  entering 
AD  and  leaving  BC  is  mixed  fluid  that  is  a  consequence 
of  a  "mixer"  in  .  While  it  is  helpful  to  think  of  the 
mixer  as  a  vortex,  this  assumption  is  not  essential  to  the 
analysis.  The  change  in  energy  flux  between  AD  and  BC 
will  be  proportional  to  the  energy  in  the  mixer,  since  the 
increase  is  due  solely  to  a  thickening  of  the  mixing  layer, 
when  the  flow  on  AD  and  BC  is  similar.  This  has  been 
found  to  be  true  in  experiments  if  one  is  in  the 
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asymptotic  mixing  regime.  If  E  is  the  energy  per  unit 
volume,  then  conservation  of  energy  in  the  moving 
frame  allows  us  to  write 

^  {  )  EdV  .  I  (E^  +  E^idV  }  = 

Vo 


M  -2ph(U  -  U  )  +  (U,  -  0  )  +  (2) 

I  c  111c 


Poh„ (U„  -  U  ) 1 

*^2  2  '  2  c'  J 


(  EdV  -  M  { 


EdV  «  M  ^-2ph(U  -  U  )  + 


M  a.  o  y 

=1  ,  7,-1  2 

|l  +  -  M 


"l  -  ^ 


2  c. 


”c,^2P2^2  ,  7,-1 

2  1  ^  12  „  2]  ^ 


2  c. 


(5) 


where 


1  2 
E  =  pe  +  2 


ph  =  p  +  E 
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Re-employing  Equation  (4)  to  relate  p.hj  to  pj  in  the  last 
(3 )  two  terms  on  the  right  hand  side,  we  write 


AA 


Here  q  is  the  velocity  relative  to  .  Throughout  we  will 
assume  that  the  dominant  part  of  the  kinetic  energy 
I  q  I  ^  is  given  simply  by  the  square  of  the  streamwise 
velocity.  This  is  exactly  true  outside  the  mixed  region 
and  approximately  so  within  the  layer,  if  the  usual 
boundary  layer  assumptions  are  made.  Thus 

q2  =  (U  -  U/  . 

Use  of  the  perfect  gas  relation  allows  us  to  write 


p^h. 


(4) 


Then  noting  that  E.| ,  Ej  are  constant,  and  dV.|  /dt  =  - 
AAU^,  Equation  (2)  becomes 
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This  expression  simplifies  greatly  if  we  invoke  the 
assumptions  explicit  to  the  convective  Mach  number 
concept,  namely  equal  static  pressures  (pj  =  Pj )  and 
equal  specific  heat  ratios  (7,  =  72  ).  This  allows  us  to 
combine  terms  and  explicitly  calculate  pR  using  the 
Crocco-Busemann  relation  (assuming  the  turbulent 
Prandtl  nurr.ber  is  1,  see  Appendix  A  for  details).  This 
yields 
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ENERGY  AND  MIXING  OF  AN  ISOLATED, 
CONVECTING  VORTEX 

Consider  the  isolated  vortex,  shown  in  Figure  4, 
convecting  in  a  uniform  free  stream  of  Mach  number 
M^.  The  mixing  rate  per  unit  span  of  such  a  two- 
dimensional  vortex  can  be  defined  approximately  as 


M  ‘Yp 

^  -fzr  I'  *  ^  -  *21 


where  v  is  the  perturbation  velocity  normal  to  the  x-axis, 
is  the  density  of  the  free  stream,  and  x^  is  chosen  to 
be  large  enough  to  encompass  the  major  influence  of  the 
vortex's  velocity  field.  The  perturbation  energy  of  this 
vortex  fe,  is  defined  as 
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Under  these  assumptions  the  increase  in  energy  of  the 
mixed  region  exactly  equals  the  decrease  in  energy  of  the 
two  streams. 

Defining  a  mixing  rate  m  as  the  rate  of  increase  of 
mass  in  the  mixed  volume  V„ 


n,  =  _  1  pdV  =  =  2pMU^  ,  (8) 


where  V  is  some  suitably  chosen  control  volume,  and  Aq 
is  the  perturbation  velocity  due  to  the  vortex,  related 
approximately  to  its  strength  r  by  T  “  2)rAqrj,  (r^.  =  core 
radius).  It  can  be  argued  strictly  on  dimensional 
grounds  that  r  is  proportional  to  the  root  of  the  energy 
Ej^  ^  ^  .  A  similar  proportionality  holds  for  the  mixing 
rate  m,  by  identifying  the  perturbation  velocity  v 
approximately  with  Aq.  As  shown  in  Appendix  B  this 
approximation  is  correct  provided  the  perturbation 
Mach  number  satisfies 
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If  we  then  define  a  nondimensional  mixing  rate  m 


we  can  integrate  Equation  (7)  to  obtain  the  growth  of  the 
mixing  region's  energy  in  terms  of  the  mixing  rate.  Thus 


m  =  m/p^U„S,  (12) 

where  S  is  the  area  through  which  fluid  is  transported, 
we  are  able  to  write  a  scaling  relation 


E  =  EdV  =  C  Tm(l  +  ^ M  )t  ,  (9) 

S  /  V  ^  c 
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where  T  is  an  average  temperature  defined  by  T  =  p/Rp. 

Note  that  this  expression  is  independent  of  the 
mechanism  for  mixing.  In  the  following  section  we 
derive  another  relation  between  mixing  rate  and  energy, 
but  this  for  a  specific  mixing  model,  the  vortex. 


Here  m^  is  the  mixing  rate  associated  with  a  reference 
vortex  of  energy  6^^ .  Equation  (13)  is  the  second  relation 
between  mixing  rate  and  energy,  but  applied  to  a  steady 
state  vortex.  Equation  (9),  on  the  other  hand,  relates  an 
increasing  energy  to  a  fixed  mixing  rate  and  an 
expanding  control  volume.  However,  because  of  the 
slow  growth  of  the  control  volume,  we  may  assume 
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Equation  03)  holds  in  the  quasi-steady  sense  and 
associate  m  (=  rn/  S)  and  there  with  the 

corresponding  quantities  in  the  energy  balance.  Equation 
(9).  This  allows  us  to  eliminate  from  Equation  (9)  and 
write 

-  2 
m 
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The  quasi-steady  assumption  is  equivalent  to  replacing 
Ej.,  Ejj  in  Equ^ion  (13)  by  their  nondimensional 
counterparts  E^,  E^  where 


-  ^o 

^o  E  p  V 

CO  O 


(It) 


E  =  CT  +  iu^  = 

o»  V  <o  2  00 


C  T  (1  + 

V  00 


7 (7-1) 
2 


(16) 


and  V  is  now  recognized  as  the  slowly  growing  control 
volume  of  the  original  energy  balance.  Equation  (2). 
This  returns  Equations  (13)  and  (14)  to  a  form  where 
mixing  rate  m  is  constant. 

What  is  the  most  appropriate  choice  for  the  scale 
quantities  T^  and  in  Equation  (16)?  In  the  frame  of 
reference  moving  with  the  vortices  the  mixed  layer  has 
temperature  T,  and  the  oncoming  free  streams  are  at  a 
Mach  number  of  .  Picking  these  as  reference  values, 
and  also  choosing  U^S  in  Equation  (12)  such  that 

UoS  =  V  ,  (17) 

we  obtain 

-  2 
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m  =  .  (18) 
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an  isolated  vortex.  In  the  next  section  we  obtain  this 
variation  analytically  after  a  novel  substitution  which 
replaces  the  reference,  isolated  vortex  with  the  bound 
circulation  of  a  thin  lifting  airfoil. 

THE  AIRFOIL  ANALOGY 

As  shown  in  Appendix  B,  the  energy  of  a 
vortex  is  proportional  to  its  rotational  kinetic  energy 
(RKE),  provided  any  temperature  perturbations  are  due 
solely  to  the  velocity  field.  Thus  we  may  rewrite  the 
equivalence  relation.  Equation  (18),  as 
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where  jj  is  some  proportionality  constant. 

If  we  now  wish  to  obtain  the  ratio  of  nondimensional 
mixing  at  to  its  incompressible  counterpart  (M^  =  0), 
we  write 


m 

tr . 

1 

where  the  subscript  i  stands  for  incompressible 
conditions.  It  is  this  quantity,  plotted  from  the 
experimental  data,  which  most  clearly  shows  the  effects 
of  decreased  mixing  as  increases. 

Lcle,^  among  others,  has  noted  that  the  flow  field 
around  the  rollers  in  a  two-dimensional  shear  layer  for 
M,.  >  .6  bears  a  strong  resemblance  to  that  about  a 
transonic  airfoil.  We  can  obtain  an  analytic 
representation  of  the  right  hand  side  of  the  above  mixing 
ratio  equation  bv  invoking  this  similarity.  Assume  that 
the  isolated  reference  vortex  in  a  freestream  of  M  =  M 

C 

can  be  replaced  by  the  bound  circulation  produced  by  a 
thin  lifting  airfoil  at  constant  angle  of  attack.  Then  in  the 
context  of  thin  airfoil  theory  we  may  solve  the  Prandtl- 
Glauert  equation  for  both  the  mixing  m^, ,  defined  similar 
to  Equation  (10),  and  the  rotational  kinetic  energy 
(RKE)jj ,  as  a  function  of  the  Mach  number.  The  details 
are  shown  in  Appendix  C.  The  result  for  the  mixing 
ratio  is 


This  equation,  one  of  the  primary  results  of  the  paper, 
shows  that  the  nondimensionalized  mixing  rate  of  the 
shear  layer  m  can  be  related  to  the  nondimensionalized 
mixing  m^  of  an  isolated  vortex,  of  energy  E^ ,  fed  by  an 
oncoming  flow  at  the  convective  Mach  number  .  This 
allows  us  to  obtain  the  variation  of  shear  layer  mixing 
with  M^  by  examining  the  corresponding  variation  for 


2 

n,.  (1  -  M  /2) 

1  c 

This  relation  is  plotted  in  Figure  5,  along  with  the  data  of 
Papamoschou  and  Roshko,  their  transcription  of 
Bogdanoff's  own  compendium,  and  two  points  from 
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Samimy  and  Elliott.  Papamoschou  and  Roshko's  earlier 
data  (AIAA  86-0162)  were  based  on  visual  thickness 
judged  from  schliercn  pictures,  while  their  later  data 
comes  from  pitot  surveys  of  the  shear  layers.  For 
convective  Mach  numbers  less  than  .8,  Equation  (21) 
captures  the  correct  qualitative  trend  of  mixing  decrease 
and  goes  right  through  the  experimental  scatter. 
Numerical  integration  of  the  airfoil  problem  using  the 
nonlinear,  transonic  small  disturbance  equation  for  the 
higher  range  of  Mach  numbers  does  not  substantially 
alter  the  shape  of  the  mixing  ratio  curve,  except  very 
near  M^  =  1 .  This  effect  can  be  seen  as  the  dotted  curve 
in  the  figure. 

What  are  the  physical  effects  embodied  in 
Equation  (21)?  One  can  immediately  identify  the 
downwash  due  to  the  airfoil  with  the  entrainment 
velocity  of  the  layer.  For  simplicity  of  analysis  we  have 
kept  downwash  constant  and  shown  that  the  rotational 
kinetic  energy  due  to  this  downwash  decreases  with  M^ . 
In  a  mixing  layer  at  a  given  velocity  ratio  it  is  perhaps 
more  physical  to  say  that  rotational  kinetic  energy  is 
kept  constant  (being  essentially  a  function  of  the  velocity 
difference)  and  that  entrainment  velocity  (downwash) 
decreases  with  M^ .  Either  way  the  decrease  is  seen  to  be 
a  function  of  the  compressibility  parameter  1  -  M^ .  The 
largo  structures  have  an  induced  velocity  field  which 
attempts  to  deflect  fluid  from  the  freestreams  into  the 
layer.  As  convective  Mach  number  increases,  the  energy 
of  this  field  goes  increasingly  into  compression  of  the 
fluid  as  it  is  deflected,  and  not  just  into  deflection,  as  it 
docs  for  the  incompressible  case. 

Our  physical  explanation  of  the  mixing  decrease  is 
at  variance  with  that  suggested  previously.^  We  agree 
that  compressibility  is  the  source  of  the  decrease 
(baroclinic  torques  play  no  part),  but  disagree  on  how  its 
effects  are  manifested.  It  is  important  to  remember  that 
the  airfoil  analogy  employs  potential  theory  to  obtain  its 
analytical  expression  for  the  mixing  decrease.  Thus  it 
ascribes  decreased  mixing  to  irrotational  effects.  This  is  in 
sharp  distinction  to  Lele's  explanation,  which  is  based  on 
arguments  concerning  vorticity  redistribution. 

THREE-DIMENSIONAL  EFFECTS: 

SWEPT  OR  BIFURCATING  VORTICES 

The  predictions  of  Equation  (21)  for  mixing  ratio 
decrease  clearly  break  down  for  M^.  >  .7-.8.  It  is  in  this 
Mach  number  range  that  two-dimensional  numerical 
simulations  show  shocks  forming  on  the  spanwisc 
vortices.  However,  experimental  evidence  for  the 
appearance  of  such  shocks  in  real  shear  layers  is  lacking. 
In  recent  work  by  Papamoschou,^  ’  which  postulated 
the  existence  of  shocks  to  explain  discrepancies  between 
measured  and  theoretical  values  of  M  ,  visual  evidence 
for  them  in  Schliercn  pictures  could  not  be  found. 

It  is  our  conjecture  that  for  M^  >  .7-. 8  three- 
dimensional  flow  becomes  increasingly  important  in  the 


mixing  layer,  and  that  the  flow  adjusts  itself  through  this 
third  dimension  so  as  to  avoid  shock  formation.  We  do 
not  know  if  this  can  be  enunciated  as  a  general  principle, 
but  we  can  point  immediately  to  a  celebrated  example, 
the  transonic  area  rule  due  to  Whitcomb,’  ^  in  which 
allowance  for  cross  flow  delays  appearance  or  weakens 
shocks.  In  the  case  of  the  shear  layer  we  suggest  two 
possible  scenarios  for  three-dimensionality:  vortex 
sweep  or  vortex  bifurcation.  In  either  case  inclusion  of 
these  effects  in  the  original  analysis  predicts  leveling  off 
of  the  mixing  ratio  curve  at  the  critical  Mach  number 
where  three-dimensionality  begins.  The  theory  does  not 
predict  what  this  Mach  number  will  be;  it  must  be  gotten 
from  experiment  or  simulation. 

Consider  the  swept  vortex  arrangement,  shown  in 
plan,  on  Figure  6.  Here  the  mixers  are  aligned  at  some 
angle,  $,  to  the  cross  stream  direction.  Figure  3  is  still 
valid  if  the  x  coordinate  is  aligned  normal  to  the 
structures  rather  than  to  the  shear  layer  flow  direction. 

Since  the  quantity  m  in  the  control  volume 
analysis  is  nondimensionalized  by  the  plan  area  S,  the 
analysis  for  the  swept  vortex  is  the  same  as  for  the 
unswept  case  except  that  the  reference  vortex  related  to 
m^  is  now  in  a  flow  with  Mach  number  M^  and  swept  at 
angle  6  to  this  flow.  Since  from  swept  wing  theory  only 
the  Mach  number  normal  to  the  vortex  is  important  in 
calculating  I  vl,  m^  is  calculated  for  a  vortex  in  a  flow 
with  a  Mach  number  of  M^cosfl.  If  M*  denotes  the 
critical  value  of  M^  at  which  shock  waves  start  to 
appear,  and  if  the  vortex  sweeps  so  as  to  avoid  shock 
waves,  then 

M'  =  M|.cose  ,  (22) 

and  it  follows  that  both  m^  and  (RKE)^  in  Equation  (20) 
will  remain  fixed  at  their  value  at  M  .  Consequently 
m/m-  will  remain  frozen  at  its  critical  value  (Mj.=:Mi 
when  Mj.  >  M*.  The  sweep  assumption  may  break 
down  for  Mach  numbers  much  greater  than  M*  because 
of  nonlinearity,  but  it  is  expected  that  the  analysis  given 
above  will  be  valid  for  low  supersonic  values  of  M^. .  The 
value  of  M  remains  to  be  decided.  If  a  value  of  0.77  is 
taken,  then  the  results  shown  in  Figure  5  as  the  chain 
dash  line  can  be  obtained.  This  value  is  in  the  range 
predicted  by  numerical  simulations  0.7  <  M  <  0.8  and  is 
the  subject  of  further  study. 

A  similar  effect  is  calculated  if  it  is  assumed  that 
the  spanwise  vortices  bifurcate  into  a  mixture  of 
spanwise  and  streamwise  vortices  such  as  sketched  in 
plan  view  on  Figure  7.  This  geometry  could  simply  be  a 
strengthening  of  the  roller-braid  structure  observed  in 
subsonic  mixing  layers.  As  is  shown  in  Appendix  D,  the 
total  energy,  E,  in  this  system  can  be  found  from 
Equations  (D1-D4)  and  gives 
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where  the  caret  (  )  refers  to  the  streamwise  vortex,  the 
tilde  (  )  to  the  spanwise  vortex,  and  o  is  a  measure  of  the 
partition  of  energies  between  the  two  vortices.  The 
velocity  in  the  spanwise  vortex  is  proportional  to  the 
strength  of  the  vortex,  which  in  turn  is  proportional  to 
£1  /  2  .  If  the  peak  velocity  in  this  vortex  is  to  be  held  at 
sonic  speed  (critical  value),  then  it  follows  that 


Consequently,  Equation  (D16)  can  be  written  as 
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where  denotes  a  value  at  the  convective  Mach 
number  M*  at  which  the  spanwise  vortex  has  a  sonic 
peak  velocity. 

From  Equation  (D4) 

aE*=E,  (25) 

and  hence  from  Equations  (23)  and  (24) 
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Using  Equation  (D4)  it  follows  from  Equation  (26)  that 


Equations  (28)  and  (30)  are  two  equations  for  m/ mj  and 
a  in  terms  of  m‘/mj.  If  the  critical  Mach  number  is 
taken  to  be  0.77,  then  this  theory  gives  the  result  shown 
in  Figure  5  as  the  dashed  line. 

CONCLUSIONS 

We  have  presented  a  theory,  free  from  linear 
instability  arguments,  which  quantitatively  predicts  the 
decreased  mixing  rates  observed  in  compressible  free 
shear  layers.  The  primary  mechanism  underlying  this 
decrease  is  the  attenuating  effect  compressibility  has  on 
the  large  scale  structures'  ability  to  deflect  fluid  away 
from  the  free  streams  and  into  the  mixing  layer.  The 
flow  energy  available  for  entrainment,  which  comes 
from  the  induced  velocity  field  of  the  structures,  is  now 
partitioned  between  deflection  and  compression  work 
that  occurs  during  deflection.  This  is  an  irrotational 
effect.  Thus,  there  is  no  connection  between  this 
mechanism  and  the  attenuation  of  vorticity  re¬ 
distribution  scheme  proposed  by  Lele.  In  incompressible 
layers,  "pushing"  on  the  free  stream  fluid  moves  it  into 
the  layer;  in  compressible  layers,  "pushing"  results  in 
some  compression  also,  and  thus  less  deflection.  This 
phenomenon  is  explicit  in  the  airfoil  analogy,  where  a 
constant  downwash  results  in  decreasing  rotational 
kinetic  energy  as  M^  increases  (or  alternatively,  if 
rotational  energy  is  held  constant  as  M^  increases,  the 
downwash  must  decrease). 


or,  using  Equation  (18) 


where  m^  is  the  value  of  m  at  incompressible  speeds. 

Since  the  spanwise  vortex  does  not  change 
character  with  Mach  number,  it  follows  that  in  Equation 
(D16)  (using  Equation  (18)) 


Reduced  linear  instability  is  not  the  cause  of 
reduced  shear  layer  mixing.  Rather  both  are 
manifestations  of  the  retarding  effects  of  compressibility. 
To  enhance  mixing  one  must  decrease  the  effects  of 
compressibility,  not  attempt  to  increase  linear  instability. 
Even  if  initial  instability  were  somehow  enhanced,  the 
long  term  growth  of  structures  would  still  be  subject  to 
the  attenuating  effects  of  compressibility.  This  is 
probably  one  of  the  reasons  why  trailing  edge  treatments 
have  little  or  no  effect  on  compressible  mixing. 

Compressibility  effects  can  be  reduced  by 
assuming  the  flow  becomes  three  dimensional.  Our 
theory  suggests  that  adjustment  of  the  flow  in  this 
manner,  in  a  way  that  avoids  shocks  at  high  values  of 
Mj.,  explains  the  levelling  of  the  growth  rate  curve  for 
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>  .7-. 8.  Indeed,  it  appears  that  Nature  has  already 
found  a  mechanism  to  keep  mixing  finite.  Without  the 
alternative  of  three-dimensional  flow,  mixing  would 
decrease  still  further,  the  energy  available  for  mixing 
having  been  further  dissipated  by  shocks.  A  difference 
in  the  amount  of  flow  three-dimensionality  from  facility 
to  facility  may  explain  some  of  the  experimental  scatter 
in  the  literature.  There  is  certainly  wide  variation  in  the 
aspect  ratio  of  the  test  sections  of  the  various  tunnels. 
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APPENDIX  A 


The  Crocco-Busemann  relation  for  turbulent  h  -  C  T  +  —  (A-l) 

Prandtl  number  Pj  =  1  regulates  the  temperature  field  in  ^ 

a  boundary  layer  to  be  strictly  a  function  of  the 

streamwise  velocity  U'  we  can  obtain  h(y).  Then  to  evaluate  the  average  value 

of  pK  across  the  layer  we  write 


(U'-c  ) 

c  T  - - -  +  C„  , 

p  2  2 


H  I 


P  (y)h  (y)dy 


where  with  boundary  condition  T  =  Tj  at  U'  =  U'j ,  and 


T  =  Tj  at  U'  =  U'j  we  obtain 


U'  +  U'  C  (T  -  T  ) 

C,  =  — ^ - -  +  -P — I - — 

1  2  (U'^  -  U'^l 


1  7P 
2L  7  -  1 


U'  -  U'  2  ' 


(Acosh  y  +  2Bsinh  y  cosh  y 

2  2  ^ 
acosh  y  +  2bsinh  y  cosh  y  +  csinh  y 


where 


C  (T,  -  T^l  ,2, 


*  "J’  *  1  1  1 


1  2' 


C  (T-  +  T  )  0  ^  +  U  ^ 

A  =  -£.-2-  .  ■>_ 


If  we  assume  a  velocity  profile  in  the  frame 


U'  +  U'  U'  -  U' 

0'  (y)  =  - - -  +  - - -  tanh  y  " 


C  (T  +  T  )  (U  ^  -  u/) 

B  =  -P_i - i_  +  - i— 

4  8 


(U'  -  U'  )  +  (U'  -  U'  J 

1  C  Z  C 


C  (T.  -  T„) 

pi  z 


(U'  -  a'  )  -  (U'  -  U'  ) 

1  c  2  c 


tanh  y  (A3) 


2  I  2 


^1  +  U^  Uf  “2 

- r -  +  - T -  tanh  y 

2  2-' 


This  integral  may  be  evaluated  with  the  aid  of  formula 
2.452.2,  page  111,  Gradshteyn  and  Ryzhik,  "Table  of 
Integrals,  Scries,  and  Products,"  Academic  Press,  1980. 


U,  =  U' ,  -  U  etc.  , 
1  1  c 


we  can  calculate  the  tem{>erature  distribution,  and  from 
that  the  density  distribution  using  the  perfect  gas  law. 


Similarly,  using  the  relation 
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ps  =  1  —ip - 1 _ 

^  4h^  -  (a  +  c)^ 


U'  -  U  2  „  U'  -  U  2 

^  ■  l-^l  ■  “/ -  l-^l  -  .-%  • 

ai  1  2  a^ 


[4Bb  -  A(a  +  c) ]2L  + 


[Ab  -  B (a  +  c) ] log 


2  c 


a  +  2btanh  L  +  ctanh  L 
2~ 

a  -  2btanh  L  +  ctanh  L 


2Ab  -  2B(a  -  c)  -  Ac (a  +  c)  , 

+  -  log 

-  72  * 

2'|b  -  ac 


Note  that  this  result  is  largely  independent  of  the 
assumed  form  of  the  velocity  profile  between  the  two 
streams,  provided  L  is  large  compared  to  the  layer 
thickness. 


APPENDIX  B 


Ictanh  L  +  b-  'lb  -  ac 


The  energy  due  to  the  vortex  is  given  by 
E  =  C,(T-Tj  +  -^(q-qj2  , 


(ctanh  L  +  b  +  Nb  -  ac 


where  denotes  ambient  conditions  and  q  is  the 
velocity. 


(-ctanh  L+b+ib  -ac 


From  the  isentropic  relations  the  temperature  due  to  the 
local  Mach  number  of  the  vortex  is 


(-ctanh  L  +  b-'lb  -  ac 


As  L  becomes  large  the  second  and  third  bracketed 
terms  go  to  finite  values  and  thus  make  negligible 
contribution  to  pfi.  Evaluation  of  the  first  term  gives 


^  |l  +  -^ 
-11  4C 


2  2 

"l  ,  "2 

C  T,  4C  T„ 
pi  p  2 


1=  [1  + 

T  2 


where  is  the  stagnation  temperature  of  the  vortex 
(i.e.,  when  M  is  zero),  and  M  is  the  Mach  number  given 
by 


.2  . 


2  2 
U  U 

7-11  7-12 


12  +  I--  I  J-  +  2_l2  _i_  I 

7  -  1  1  4  2  4  2  J 


However, 


Aq  =(q-qj  , 


a^  =7RT. 


Equation  (B2)  can  be  expanded  to  give 
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APPENDIX  C 


(B6) 


so  to  a  first  approximation 


T  _  ,  _  (7-1)  Aq 

T  2  2' 

00  a 


where  a  is  the  local  speed  of  sound.  Also 


(B7) 


2  T 

^  — 


7RT  [1  -  M^- 


a  2  [1  -  „2. 


1 

(B8) 


To  a  first  approximation  in  Aq^  then 
1I_  -  -  (7-1)  Aq^ 

(t  J  22' 

00  a 

00 

Equation  (Bl)  can  now  be  written  as 


T  112 

E  =  c  T  ^  -  1  +  ^  Aq 

V  00  1  12 


=  c  T  f-  ^ J  .  I  Aq^ 

V  00 1  2  2  J  2 

a 


(B9) 


(BIO) 


Aq  (7-1) 
2  7 


Analysis  for  Subsonic  Flow 
Consider  the  Prandtl-Glauert  Equation 


^6  +6  =0  , 
^xx  ^yy  ' 


(Cl) 


where  ^  is  a  perturbation  velocity  potential,  such  that, 

u  =  (C2) 

v  =  ^' 

and  u  and  v  are  normalized  perturbation  velocities  in  the 
X  and  y  directions  respectively,  x  and  y  are 
nondimensionalized  to  the  airfoil  chord,  2x^  ;  p  is  given 
by 

^  =  (1-M/)1/2  (C3) 

Equation  (Cl)  can  be  transformed  by  the  substitution 

y  =  ^y  (C4) 

to  give  Laplace's  equation 


6  +6 
^xx  ^yy 


:0. 


(C5) 


For  the  flow  over  an  airfoil  with  lift  the  classic  thin  airfoil 
analysis  gives 


Au  (x 


-'■r  j 


Au (X,  o)  y 

2  -2 
(x-n  +  y 


(C6) 


i  ( 

2»  jo 


Au  (x,o)  (x-t) 

2  -2 
(x-{)  +  y 


d« 


(C7) 


where  Au  and  are  defined  by 


Au(x,y)  =  J-  [u(x,y)  -  u(x,-y  )1 


(C8a) 


Hence  E  is  proportional  to  Aq^  .  Since  Aq  is  a 
perturbation  velocity  it  is  proportional  to  the  strength  of 
the  vortex,  P,  and  hence 

EecI^ 

Since  the  rotational  kinetic  energy  is  proportional  to  Aq  ^ 
it  follows  that  E  is  proportional  to  the  rotational  kinetic 
energy. 


=  r  lv(x,y)  +  v(x,y)l . 


For  a  given  angle  of  attack, 


or 


V  =  -Oq  =  , 


Py=-ao/P  ■ 


(C8b) 


(C9) 


(CIO) 


By  putting  y  =  0  in  Equation  (C7)  it  may  be  seen  that 


-n  - 


Au(x,o)  =  AUj(x,o)/^  , 


(Cll) 


whi'ro  the  subscript  i  denotes  a  value  at  incompressible 
speeds,  that  is,  at  p  =  1.  It  may  also  be  seen  that 

v(x,y)  =  =  V  j(x,y )  (Cl  2) 

and  that,  from  Equation  (C6), 


1  (  2  -  - 

-  p  Au  (x,y)dxdy 

2  00  J  s  X 

1  (  -  2 

-  p  V.  (x,y)dxdy 

2  CO  J  s  1 


(C19) 


Au(x,y)  =  ^Uj(x,y) .  (C13) 

If  the  mixing,  m,  is  defined  as  the  movement  of  fluid 
from  the  upper  half  plane  to  the  center  and  vice-versa, 
then 


which  are  independent  of  Mach  number. 

Using  Equations  (Cl  5)  and  (Cl  8)  in  Equation  (20) 
of  the  main  text  gives 


m  =  p 


1 

Ivl  „  dx  +  Ivl  dxl  ,  (C14) 

1  y=0  J  y=0  J 

■^1  1 


which,  using  Equation  (Cl  2)  gives 
m  =  mj , 


where  mj  is  the  mixing  for  incompressible  flow  for  the 


The  rotational  kinetic  energy  (RKE)  is  defined  in  a 
small  disturbance  approximation  by 


RKE  “  4  p  +  V  ^1  dxdy  ,  (C16) 

2  ee  S 


m 
m . 


h  *  *2 

■l  ♦ 


(C20) 


If  the  airfoil  is  small  compared  to  the  extent  of  domains, 
_  _  there  is  no  difference  between  Au(x,y)  and  V(x,y),  that  is, 
’  the  airfoil  resembles  a  point  vortex.  In  this  case 


and  Equation  (C20)  becomes 


m. 

1 


1  -  M  /2 
c 


(C21) 


(C22) 


where  S  is  half  of  the  flow  domain. 

Equation  (C16)  can  be  written  as 


APPENDIX  D 


Effect  of  Multiple  Vortices 


1  ^at  2-2 

RKE  =  J  —  (x,y) /P 


(CIV) 


by 


From  Equations  (9)  and  (13)  the  energy  E^  is  given 


+  (x,  y)  ]  dxdy 


^  —  —  r  2 

E  =  mC  T  1  +  ^  „ - -  M  p  U  St  (Dl) 

s  vl  2  cjoo 


RKE  - 


1 


where  in  is  the  mixing  due  to  the  structures.  Let  the 
mixing  be  caused  by  N  vortices  on  the  mixing  line,  each 
(Cl  8)  with  a  mixing  rate  of  mj  such  that 


Here 


(D2) 
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Using  the  scaling  of  Equation  (13)  we  obtain 


Again,  using  the  scaling  of  Equation  (13) 


where  the  subscript  "o"  denotes  a  mixer  with  energy  E^ , 
and  Ej  j  is  the  energy  associated  with  the  jth  structure.  If 
the  energy  of  each  vortex  is  a  fraction  of  the  total 
energy  E^,  then 


(D9) 


and 


(D4) 


It  is  assumed  that  the  mixing  is  proportional  to  the 
square  root  of  the  energy  of  the  mixer  as  is  the  case  for 
the  vortex. 


N 

^  a.  -  1  .  (D5) 

j-1  ’ 


Using  Equation  (D2)  -  (D4)  in  Equation  (Dl)  gives 


Now  assume  that  the  energy  of  each  structure  is  a 
fraction  aj  of  the  total  energy  E^  and  Equations  (D4)  and 
(D5)  are  satisfied.  Use  Equations  (D8),  (D9),  and  (D4)  to 
eliminate  Eg.  Thus 


m 


m  a  . 
o  3 

E 

o 


(06) 


where 


E  =  E  / 
o  o 


C  Tp  U  St  1  + 
V  o  o 


Thei4ach  number  dependence  only  arises  in  the 
term  m^  /E^  ^  ^  ^  and  hence  the  mixing  will  decay  with 
Mach  number  no  matter  how  many  vortices  there  are. 
However,  as  before,  this  analysis  neglects  the  viscous 
core  of  the  vortices.  If  the  vortices  arc  small,  the  viscous 
core  is  larger  relative  to  the  inviscid  part  than  in  larger 
vortices,  so  overall,  the  mixing  may  not  decrease  as 
much  with  Mach  number. 

Consider  now  the  case  where  there  arc  two  types 
of  mixing  structures,  with  mixing  denoted  by  m  and  m 
such  that 


N  m  T 

£  »~r/2  “j  H  “ 


j  ^  1/2 
j=N+l  ■'  E 


,  (DIO) 


where 


E  =E/|cTpUSt|l+ 
o  o  I  V  o  o  ' 


7(7-1)  „  2 


“/ll 


(Dll) 


E  =E/|cTpUSt[l  +  M  ^  1  I  . 

o  o  I  V  '^o  o  I  2  c  J  J 

In  order  to  alter  the  effect  of  Mach  number  on  the 
mixing,  the  typical  term 


E 

o 


N 

m  =  'y  m .  -t 
j-1  ^ 


T  . 

Z 

j-N+l-* 


must  be  less  dependent  on  the  Mach  number  than  is  the 
case  for  the  simple  vortex. 

( D8 )  Note  that  the  structure  denoted  by  the  caret  can  be 

regarded  as  the  viscous  core  and  that  denoted  by  a  tilde 
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can  be  regarded  as  the  inviscid  part  of  the  vortex  which 
decreases  rapidly  as  M_,  — 1. 

If  the  structure  denoted  by  a  caret  is  a  vortex  with 
a  streamwise  axis,  then  the  transonic  small  disturbance 
equation  gives  the  potential  associated  with  the  vortex  as 

=0.  (D12) 


where 


k  =  [3  +  (7-2)M/]M^2  (Q13) 

If  the  vortex  is  uniform  or  nearly  uniform  along  its 
axis,  then  the  x  derivatives  in  Equation  (D12)  are 
negligible  and 


^yy  ^zz  ~  ^  '  (D14) 

which  has  no  dependence  on  M^..  Thus  the  mixing 
associated  with  a  vortex  aligned  in  the  streamwise 
direction  does  not  decay  as  — >  1.  Furthermore, 

there  are  no  shock  waves  associated  with  this  vortex  and 
hence  no  increase  in  entropy. 

In  contrast,  for  a  spanwise  vortex  the  potential 
equation  is 


(1-M,2  +«yy  =0,  (D15) 

which  can  give  shock  waves  and  hence  an  entropy 
increase. 

Consider  now  the  mixing  given  by  Equation  (DIO) 
for  a  single  vortex  and  a  single  streamwise  vortex.  Thus 


-  1/2 
E 


a  +  (1  -  Of) 


1/2 


1/2 


(016) 


where  a  is  a  measure  of  the  partition  of  the  energy 
between  the  vortices. 

It  is  a  reasonable  hypothesis  that  the  flow  will  try 
to  develop  so  that  entropy  increases  are  kept  to  a 
minimum.  Hence,  as  increases,  the  fraction  a  will 
get  smaller  and  the  size  of  the  spanwise  vortex  will 
decrease  while  the  streamwise  vortex  will  increase.  This 
effect  is  in  addition  to  the  inviscid  decay  of  the  spanwise 
vortex  noted  earlier. 


-14- 


Figure  2.  Initizd  spreading  rates  of  forced  FIGURE  3.  Control  volume  for  energy  bal 
two-dimensional  shear  layers  at  various  con-  ance  in  a  shear  layer, 
vective  Mach  numbers 


Figure  4.  Isolated  vortex  model 


Figure  6.  Plan  view  of  swept  vortices 


Figure  5.  Experimentally  measured  layer 
growth  rates  compared  to  two-  and  three- 
dimensional  theory. 


Figure  7.  Plan  view  of  mixture  of  spanwise 
and  streamwise  vortices 
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